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<N ■ ABSTRACT 

O 
D 

Q ■ Magnetic reconnection is a fundamental process in space and astrophysical 

plasmas in which oppositely directed magnetic fields changes its connectivity 



and eventually converts its energy into kinetic and thermal energy of the plasma. 
Recently, ubiquitous jets (for example, chromospheric anemone jets, penumbral 
microjets, umbral light bridge jets) have been observed by Solar Optical Tele- 
scope on board the satellite Hinode. These tiny and frequently occurring jets 



P^. are considered to be a possible evidence of small-scale ubiquitous reconnection in 

Q . the solar atmosphere. However, the details of three dimensional magnetic con- 

figuration are still not very clear. Here we propose a new model based on three 
. dimensional simulations of magnetic reconnection using a typical current sheet 

magnetic configuration with a strong guide field. The most interesting feature 
is that the jets produced by the reconnection eventually move along the guide 
field lines. This model provides a fresh understanding of newly discovered ubiq- 
i/-^ . uitous jets and moreover a new observational basis for the theory of astrophysical 



magnetic reconnection. 

Subject headings: Magnetohydrodynamics (MHD) — Magnetic reconnection — Methods: 
numerical 



c5 ■ 1. INTRODUCTION 

Recent observations have found ubiquitous plasma jets over the solar atmosphere in var- 
ious plasma parameters and magnetic field configuration. The solar chromospheric anemone 
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jets show a cusp- or inverted Y-shaped structure which are beheved to be a resu lt of magnetic 
reconnection between a magnetic bipole and a preexisting uniform vertical field ( IShibata et al. 
20071 ). The penumbral microjets show the apparent motion almost along the vertical guide 
field comp onents in the interlock ing-comb structure of magnetic field lines in the sunspot 
penumbra ( iKatsukawa et al. The umbral light bridge jets are ejected along the vertical 

field lines emanating from the light bridge in the sunspot umbra ( IShimizu et al.ll2009[ ). There 
are many other chromospheric jets whose footpoints are not well resolved. I t has sometimes 



been proposed that sp icules may be produced by magnetic reconnection (ISuematsu et al. 



20081 : llsobe et al.ll2008[ ). The three dimensional magnetic field configuration at the footpoints 
of these jets are still puzzling. Especially, how these jets are accelerated is a fundamental 
question in solar physics, which would also give a hint to the understanding of the origin of 
astrophysical jets. The three dimensional (3D) numerical simulation may help us to under- 
stand this. For 3D magnetic reconnection, many simulations show the r econnecting process is 
much more complicated and difficult than two dirnensional (2D) case (lYokovama fc Shibata 



19951 : IChen et al.lll999l I2OOII : lYokovama fc Shibatall200ll : llsobe et al.ll2005l : 
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Some of the 3D simulations have sh own generation of flows paral lel to magnetic fleld lines 
as a result of non-null reconnection (jPontin et al.ll2005t IUgaill2010l ). However, these parallel 
flows have not been analyzed in detail. In this paper, we analyzed these flows in detail for the 
flrst time, because these flows are important as the origin of chromospheric jets and found 
that these jets ejected from the diffusion region move along the magnetic guide fleld and its 
3D structure is similar to a fan-shape (hereafter, referred to as fan-shaped jets) which differ 
from the classical reconnection theory (generally speaking, the rec onnection je t s move along 
the ambient rnagnetic fleld, hereafter , ordinary reconnection jets) ( ISweetlll958l : |Parkerlll957l : 
Petscheklll964l : iPriest fc Forbes! I2OOOI ). In this report, we give a description and analysis for 
these results. 



2. NUMERICAL METHOD 

We perform three dimensional magnetohydrodynamic (MHD) simulation of magnetic 
reconnection using a typical force free magnetic fleld. The mathematical form of this fleld is 



5, = 



(1) 



5„ = I Bini sm{[ex - (vr - ^) A/i^]/2A/i^,} 
Bin, sin [(2^-7r)/2] 



for X < —Ahx 
for |x| < Ah^ 
for X > Ah^ 



(2) 
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= { Bini co^{[ex - (tt - ^^)A/iJ/2A/i,.} 
Bini cos [(20-7r)/2] 



for X < —Ahrc 
for \x\ < Ahx 
for X > Ah^ 



(3) 



where Bini is the initial amphtude of magnetic field. The Ah^j. is the half width of the 
current sheet. 6 is the reconnection angle shown in Figure [H and it can be four values in 
our simulation cases, i.e., vr, 37r/4, n/2 and tt/S. 

The computational box is resolved by 240 x 480 x 480 grid points, as shown in Figure [TJ 
Since the magnetic field in the current sheet is a function of x, we use a non-uniform grid 
point distribution in the x direction. The eight independent variables are density (p), veloc- 
ity {vx, Vy, Vz), magnetic field {B^, By, B^) and gas pressure (p). We studied four cases for 
different initial reconnection angles (vr/3, 7r/2, 37r/4 and vr) which are define d in Figure [H 
Using the free bounda r y condition and th e CIP-MOCCT numerical scheme (lYabe fc Aoki 



1991 



Yabe et al.l Il99ll : iKudoh et al.l Il999l ) , we solved three dimensional compressible resis- 



tive MHD equations. Note that there is no thermal conduction and gravity being considered 
in our simulations because what we are really interested in is the process of magnetic recon- 
nection. 



3. RESULTS 

Figure [2] shows the 3D visualization of the gas pressure distribution at the time = 8.5 
and time = 11.5 of our typical case (the reconnection angle is vr). We added a cross section 
at z = and a velocity field plane at x = in the lower panels, in which we can see 
the X-shaped structure of ordinary reconnection outflow on the x — y plane and the velocity 
distribution of the fan-shaped jets on the y — z plane. The plasma is heated and the magnetic 
field lines disconnect and reconnect at central diffusion region. Then we have two kinds of 
jets. The one is the ordinary reconnection outflow which is shown by the X-shaped structure 
on the X — y plane. Besides, it can be seen that there are two fan-shaped jets moving along 
the z and -z direction. The biggest feature of these jets is that they can move along the 
magnetic guide field. Since the magnetic field is shearing in the current sheet, the jets can 
move in different direction (the shearing magnetic field means the direction of the magnetic 
field lines is different if they have a different position x in the diffusion region). The shape 
of ordinary reconnection jets on x — y plane is very similar to the results in 2D reconnection 
case and the velocity is about the local Alfven speed at the time = 11.5 (the speed of the 
fan-shape jets is about half of Alfven speed). 

Figure [3] shows a cartoon for explaining the fan-shaped jets (red arrows) and the or- 
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dinary reconnection jets (blue arrows). The fan-shaped jets are ejected from the central 
diffusion region in different direction which is due to the shearing magnetic field in the cur- 
rent sheet. Finally, we have a structure which looks like a fan-shape. The reconnection angle 
for Figure [3h and b is vr, and that for Figure |3t and d is 7r/3. As we described above, the 
fan-shaped jets can move along the guide field in our typical case. For the small reconnection 
angle case, similar to the typical case, two oblique magnetic field can produce the fan-shaped 
jets in a narrow angle. This situation is very general in the solar atmosphere. If a vertical 
magnetic flux tube is twisted or sheared or if two magnetic field lines or tubes with different 
directions get close to each other (as shown in the Figure [3h and d considering x — y plane 
as the solar surface), large electric currents are generated between them, resulting in fast 
magnetic reconnection with larger amount of energy released. The ordinary reconnection 
jets and fan-shaped jets are ejected from the diffusion region simultaneously. Since the or- 
dinary reconnection jets almost parallel to the solar surface, they can only move a short 
distance and disappear. However, the fan-shaped jets can move upwards and become the 
jets along the vertical field lines. It should be noted that the density stratification effect in 
the solar atmosphere can greatly increase velocity amplitude of the slow mode shock ahead 
of the upward fan-shaped jets when it propagates to the low density region (like the upper 
chrom osphere and solar co rona) and eventually become a faster jet accelerated by the slow 



shock fIShibata et al.lll982[ ) 



In order to understand why the jets move along the guide field, we investigate the driving 
force of the fan-shaped jet using Lagrangian fluid elements (test particles). These elements 
can move with the plasma but has no effect to the result of simulations, which can show us 
the detailed information at the position of these elements in the computational box. One 
of these fluid elements is located at the edge of the diffusion region at the initial time and 
finally the element is driven as a part of the fan-shaped outflow. Its forces and velocity is 
shown by Figure S^. There are two stages for acceleration. In the first stage {time = 8 — 10), 
the Lorentz force (including magnetic tension force and magnetic pressure gradient force) 
drives this element. After that it is dominated by the gas pressure gradient in the later 
stage {time = 10 — 11.5). Hence, we conclude that the fan-shaped jets are accelerated by 
both gas pressure gradient force and the Lorentz force. That is different from the ordinary 
reconnection outflow which is only accelerated by the Lorentz force. 

Figure Hlb and c shows the dependence of velocity on the resistivity and the reconnection 
angle. All the velocity is measured at the time when reconnection rate reaches to the 
maximum for different cases. As shown in Figures |l]b, velocity are almost constant with 
different resistivity value. That is to say, there is no remarkable dependence of the velocity 
on the resistivity value. The ratio of the velocity between two kinds of jets is approximate 
0.5 (the velocity of the ordinary reconnection jets is comparable to the local Alfven speed 
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determined by the reconnecting component), which indicates that the mechanism for forming 
these two jets is different as we described above. In Figure Hh, the velocity is sensitive to 
the reconnection angle. The large reconnection angle means a strong shearing magnetic 
field lines in the current sheet. If we use the reconnection angle = as an extreme case, 
there will be no reconnection process and no jets come out. Furthermore, we found that the 
maximum speed of fan-shaped jets is still around half of the ordinary ones. 



4. DISCUSSION 

Let us briefly discuss the application of our results to various jets in the chromosphere 
as shown in Table [H Although the footpoints of these jets are not necessarily well resolved, 
we assume that the 3D reconnection occurs in the photosphere or in the low chromosphere 
at the footpoints of these jets, and fan-shaped jets are ejected from the reconnection region 
along magnetic field lines in these layers. The Alfven speed (va) in the photosphere and 
low chromosphere is about 10 km s~^ in a typical isolated flux tube outside sunspots, and is 
about 10-50 km s~^ in sunspot umbra and penumbra. Hence the va in Table [1] shows such 
local Alfven speed in the (hypothetical) reconnection region at the footpoint of these jets. 
Here, the va± is the Alfven speed based on the reconnecting component of magnetic field 
(we assume that va± = 0.2va)- The velocity of fan-shaped jets is only half of that of the 
ordinary reconnection flow (f/an ~ VA^f^). Once the fan-shaped jets are ejected, the slow 
mode MHD shock is formed ahead of the jets and propagate along the vertical magnetic 
field lines. Since the density decreases with height, the velocity amplitude at the slow mode 
shock increases with height, namely Vfan,max ~ Vjane^'^^^^ if the slow mode wave energ y 



is conserved or Vfan,max ~ Vjans'^''^^^^^ if the shock is strong (IShibata fc Suematsul Il982l ). 
where the value z is the height of the jets measured from the reconnection region and H is 
the pressure scale height. If z/H = 13.3 (assuming the scale height H is 150 km and the 
slow mode shock propagates over a height of 2000 km), we get Vfan,max ~ Vfane^'"^^^^^ ~ 
VfanG^'^ ~ 20vfan- Similar results are obtained also for the case z = 1000 km when the wave 
energy is conserved {vfan,max ~ VfanS^'^^^^)- Table [T] shows that the resulting velocity at the 
shock front {vfan,max) whcu it reaches at the top of the chromosphere is comparable to the 
actual observed velocity of these chromospheric jets. Hence, our finding of fan shaped jets 
parallel to field lines are important for understanding the origin of chromospheric jets and 
seems to be successfully applicable to ubiquitous chromospheric jets. 

We described fan-shaped jets by simulating the 3D reconnection process using a simple 
initial shearing magnetic configuration in this report. We found that the fan-shaped jets 
which are accelerated by both gas pressure gradient and Lorentz force can move along the 
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Table 1: Comparison between observation and fan-shaped jets (unit: km s 



Jets 


Observational velocity 


Va 


va± 


Vfan 




in, max 


Chromospheric anemone jets^ 


~ 10 


10 


2 


1 




20 


Spicules^ 


~ 25 


10 


2 


1 




20 


Penumbral jets^ 


50 - 150 


10-50 


2-10 


1-5 


20 


- 100 


Umbral light bright jets'^ 


28 - 180 


10-50 


2-10 


1-5 


20 


- 100 



^ Shibata et al. ( 


2007) 


' Suematsu et al 


( 


200S|) 


' Katsukawa et al. ( 


2007) 


^ Shimizu et al. ( 


2009) 



magnetic guide field lines and the velocity of these jets is about half of the local Alfven speed 
determined by the reconnecting component of magnetic field. This new finding provides us 
a new way to understanding the magnetic reconnection in 3D geometry and it is also a new 
model for explaining the solar ubiquitous chromospheric or more general astrophysical jets. 
The details of them will be studied in our future papers. 
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Fig. 1. — Initial magnetic field configuration. The box size is —4 < x < 4, —16 < y < 16, 
—16 < z < 16. The magnetic field lines rotate from x = —0.4 to x = 0.4. The 
initial uniform distributions of density and pressure are used in our simulations. The 
magnetic hydrostatic equilibrium is satisfied at the initial condition. The normalization 
units (the values used in our analysis and discussion is in nondimensional for simplic- 
ity) of length, temperature, density, pressure, time, velocity and magnetic field are Lq, 
To, Po, poToKB/m, Lq/vq, {po/ PoY^"^ and {poY^'^, respectively (kb and m mean Boltzmann 
constant and molecular mass). We assume that an anomalous resistivity, with the form 
V = i]iniCos{x7i/2Ah) cos{y'n-/2Ah) cos{zTT/2Ah) (the rjini is the amplitude of resistivity in 
the diffusion region), is localized in a small region |a;| < Ah, \y\ < Ah, \z\ < Ah, where 
Ah = 0.4. The typical magnetic Reynolds number, plasma beta and Alfven speed is 126, 
0.8 and 1.58, respectively, where the Reynolds number is calculated using rjini (0.05), Alfven 
speed (1.58) and the half length of the computational box in x direction (4). 
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Time = 8.50 b Time = 11.5 




Time = 11.5 d Time = 11.5 




Fig. 2. — Typical case of a fan-shaped jet. The 3D visuahzation (panels a-d) of the gas 
pressure distribution with the reconnection angle vr. Solid tubes are magnetic field lines and 
arrows are velocity, (a) 3D visualization of pressure distribution at the time 8.5. (b) 3D 
visualization of pressure distribution at the time 11.5. (c and d) are the same distribution 
at the time 11.5, but we added a cross section at z = (the gas pressure rendering above 
this plane is hidden) and a velocity field plane at x = 0. Moreover, (d) shows image seen 
from the right side of (c) (the view angle is different). Since the current sheet is very thin, 
a different scale in x direction is used in these panels. 
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Fig. 3. — A cartoon for explaining the fan-shaped jets, which shows both ordinary reconnec- 
tion jets and the fan-shaped jets. The blue jets are the ordinary reconnection jets and the 
red ones are the fan-shaped jets, (a) The initial magnetic configuration with the reconnec- 
tion angle tt. (b) The fan-shaped and ordinary jets with the reconnection angle vr. (c and 
d) are the same as (a) and (b) but the reconnection angle is 7r/3. Note that the magnetic 
configuration is shearing in the current sheet and the jets always move along the magnetic 
field lines, so the structure looks like a fan-shape. 
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Fig. 4. — (a) This panel shows the forces per unit mass and the velocity of one Lagrangian 
fluid element in the z-direction for this typical case. The other panels show the dependence of 
the jet velocity on the resistivity and the reconnection angle. The value Vmax is the maximum 
velocity in the computational box which indicates the speed of ordinary reconnection jets. 
vpara_max IS the maximum velocity along the initial magnetic field line which means the 
velocity of fan-shaped jets, (b) This is the dependence of the jet velocity on the resistivity 
where the rjini is the amplitude of resistivity in the diffusion region. For this case we choose 
three resistivity value 0.025, 0.05 and 0.075 with the plasma beta 1.0. (c) This panel shows 
the dependence of the jet velocity on the reconnection angle. For this case we choose four 
reconnection angle value 7r/3, 7r/2, 37r/4 and vr with the plasma beta 0.8. 



